Cosmological model with interactions in the dark sector 
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A cosmological model is proposed for the current Universe consisted of non-interacting baryonic 
matter and interacting dark components. The dark energy and dark matter are coupled through 
their effective barotropic indexes, which are considered as functions of the ratio between their energy 
densities. It is investigated two cases where the ratio is asymptotically stable and their parameters 
are adjusted by considering best fits to Hubble function data. It is shown that the deceleration 
parameter, the densities parameters, and the luminosity distance have the correct behavior which 
is expected for a viable present scenario of the Universe. 
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I. INTRODUCTION 

Recent observations of Supernova la (SNIa) suggest 
that the Universe has entered into a stage of an ac- 
celerated expansion with a redshift z ^ 1, [l|, [2|, 3- 
This has been confirmed by precise measurements of the 
spectrum of the Cosmic Microvifave Background (CMB) 
anisotropies [j|, [5[ as weU as the baryon acoustic oscil- 
lations (BAG) in the Sloan Digital Sky Survey (SDSS) 
luminous galaxy sample [6|. All usual types of matter 
with positive pressure generate attractive forces and de- 
celerate the expansion of the universe. For that reason, a 
dark energy component with negative pressure was sug- 
gested to account for a fluid that drives the current ac- 
celerated expansion (see [3| for reviews). The simplest 
explanation of dark energy is provided by a cosmological 
constant, but the scenario is plagued by a severe fine tun- 
ing problem associated with its energy scale. The vacuum 
energy density falls below the value predicted by any sen- 
sible quantum field theory by many orders of magnitude 
3, and it unavoidably leads to the coincidence problem, 
i.e., "Why are the vacuum and matter energy densities of 
precisely the same order today?" [0] . More sophisticated 
models replace the cosmological constant by a dynami- 
cal dark energy that may be a scalar field (quintessence), 
tachyon field, phantom field or exotic equations of state. 
These models fit the observational data but it is doubt- 
ful that they can solve the coincidence problem pij|.pd|. 
Besides, in a viable dark energy scenario we require that 
the energy density of the dark fiuid remains subdominant 
during radiation and matter dominant eras and that it 
becomes important only at late times to account for the 
current acceleration of the Universe. 

In several papers, it has been proposed that dark 



matter and dark energy are coupled and do not evolve 
separately [141, [IjI- The coupling between matter and 
quintessence is either motivated by high energy particle 
physics considerations [13[ or is constructed by requiring 
the final matter to dark energy ratio to be stable against 
perturbations [ll|,[lJ|. In the same sense, we will use 
two arbitrary interacting fiuids in the dark sector, with 
constant barotropic indexes, and will show that there is 
a stable solution for the ratio dark matter-dark energy, 
which satisfies the expected behavior from the matter 
dominated era until today. We explain in section II the 
general features of the model for two different interac- 
tions and find the stability conditions of the ratio dark 
matter-dark energy. The section III is devoted to obtain 
the values of parameters that best fit the experimental 
data using the recently published Hubble function H{z) 
data [l^l, extracted from differential ages of passively 
evolving galaxies. The allowed regions of probability la 
and 2(7 for the parameters are showed. Also, deceleration 
parameter, densities parameters, ratio dark matter-dark 
energy, effective equations of state and luminosity dis- 
tance are exhibited in section IV for the best fit param- 
eters, showing the correct behavior expected for a viable 
scenario. Finally, in section V we present our conclusions. 



II. INTERACTING DARK MODEL 

Let us consider a Universe modeled by a mixture of 
three constituents, namely, baryons, dark matter and 
dark energy. The Friedmann and energy conservation 
equations read 

3H^ = pi + P2 + P3, (1) 

Pl+P2+iH{pi+pi+p2+P2)^0, (2) 

P3 + 3H{p3+P3)^0. (3) 
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respectively. The subindexes 1 and 2 refer to dark matter 
and dark energy, respectively, whereas the subindex 3 to 
baryons. The baryonic matter is supposed to be a non- 



interacting component so that it decouples from equa- 
tion ^. Moreover, by considering the baryons as pres- 
sureless (ps = 0), the integration of equation Q gives 
Ps — P3 i'^0/0.) I where pg and uq represent the present 
values of the energy density and cosmic scale factor, re- 
spectively. 

Hence, equation ([2|) express the interaction between 
the dark matter and dark energy allowing the mutual 
exchange of energy and momentum. Consequently, there 
will be no local energy-momentum conservation for the 
fluids separately. However, we can decoupled equation 
([2|) into two " effective equations of conservation" as fol- 
lows 



Pi + 3H-ftpi = 0, 
P2 + iH^lP2 = 0. 



(4) 
(5) 



In the above equations it was introduced effective 
barotropic indexes 7^ (i = 1, 2) given by 



r itlpi 
Pi 



12=1-2+ lir + 



3Hp2' 



(6) 

(7) 



where r = Pi/p2 is the ratio between the energy densi- 
ties and li (i = 1,2) represent barotropic indexes of the 
equations of state pi = (7,; — l)pi. 
In addition we get the relationship 



stable whether {dAj'^/dr)r=r^ > 0. Expressing this con- 
dition in terms of Eq. (|10p . we obtain the stability con- 
dition 



(1 



-)l'f 



- F{rs) > 0, 



(13) 



where we have assumed that the barotropic indexes li of 
the two fluid are constants. Note that the simplest choice, 
a negative constant F, satisfies the latter condition. This 
fact, induce us to investigate the F constant case. So, in 
the next sections we shall analyze two different choices 
for the function F{r), namely, one refers to a constant 
and another being a variable. 



A. F{r)= constant< 

If we choose the barotropic indexes 71 and 72 as con- 
stants and the function F{r) as 



F(r) 



l + r-c 



-A7, 



(14) 



where r^ is a constant value of the ratio between the 
energy densities at infinity, then it is possible to inte- 
grate the equation ^ for r along with the effective con- 
servation equations ^ and ([5]) for the energy densities, 
yielding 



(7r-7i)^+(72'-72)=0, (8) 

and a dynamical equation for the ratio r, namely 

f = -SHrAj", AY = 7i - 72- (9) 

We assume that the effective barotropic index of the 
dark energy is given by 72 = 72 — ^('')j where F{r) is a 
function which depends on the ratio r, so that we get 



Aj" = A7 + F{r) 1 



1 



A7 = 7i-72. (10) 



By taking into account the above representation for 7I, 
Eqs. ([4]) and ^ can be rewritten as 



pi + 3i/pi7i = -'3Hp2F{r), 
P2 + 3Hp2i2 = ■iHp2F{r). 



(11) 
(12) 



Now it is possible to identify the right-hand side of each 
above equations as the energy transfer between dark mat- 
ter and dark energy. 

Now we make a detailed study of the general dynamics 
of the density ratio r as given by Eq. ^ along with Eq. 
pO)) . To this end we assume that the constant solutions 
r = rg, represents a stationary stage of the Universe. 
Then, according to Eq. ^. it means that A7'^(rs) = 0. 
Then, the constant solutions r,, will be asymptotically 



Pr=P^|.oo + (^„-^.)(^)^"[(^)^ (15) 

(16) 



/ao^'^ 



P2= P2\ 

V a 



where ag and tq are the present values of the cosmic scale 
factor of the ratio of the energy densities, whereas 

A7 r, o''°o7l+72 ,.„. 

" = -;— , /^ = 3— — . 17) 

1 + Too l+Too 

As usual we introduce the red-shift z = 1/ a— 1 and the 
density parameters at the present value zq = 0, namely 

n?=rof7°, f^^ = -4- "s^tJi. (18) 



iur 



3Hr 



with ilj + 5I2 + ^3 == 1- Once the explicit dependence of 
the energy densities are known in terms of the red-shift 
we can rewrite the Friedmann equation ([T]), thanks to 
(fTSl) and (HSl), as 



H^ = Hi {nl{l + r^)(l + zf + nl{l + zf 

+[\-nl-{\ + r^)nl]{i + zf'''), (19) 

which gives the dependence of the Hubble parameter H 
with respect to the red-shift. 



Variable F{r) 



One possible choice for the function F{r) is 



F(r}-- 



F{r) = - 






A7, 



(20) 



where Too has the same meaning as in the previous case. 
Note that the above chosen function reduces to P^ in 
the Hmit r — r^c ■ With this restriction the latter function 
satisfies the stability condition p^ . So that, the station- 
ary solution of Eq. ([5]), corresponding to the function 
(|20|) . is asymptotically stable. 

In this case the solution of Eqs. ^ and ([5]) are given 
by 

Pi-p'2^rl^ + irl~rl^){^y 



(21) 
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(22) 



(23) 



Furthermore, the Hubble parameter H can be expressed 
as functions of the red-shift, yielding 



^0 



(1 ~r/roo)(l + ro/roo) 
(1 - ro/roo)(l -t-r/roo) 



"L + (^0 



r^)(l + z)'^ +f}^(l + z)3. (24) 



Here it is interesting to call attention that the condi- 
tion Too = leads to a mixture of noninteracting fluids 
with constant barotropic indexes for both cases analyzed 
above. 



III. COSMOLOGICAL CONSTRAINTS 

The aim of this section is to determine the param- 
eters (roo,/3) for the constant F(r) case and (^00,72) 
for the variable F{r) one. These parameters, denoted 
gencrically by a^, are established through the recently 
published Hubble parameter H{z) data [l^l, extracted 
from differential ages of passively evolving galaxies. This 
H(z) function is a very interesting one, because in con- 
trast to standard candle luminosity distances or stan- 
dard ruler angular diameter distances, it is not integrated 







FIG. 1; Constant F: (3 vs. r^ 



over. The Hubble parameter depends on the differen- 
tial age of the Universe as a function of z in the form 
H{z) = —(1 + z)~^dz/dt, and it can be measured di- 
rectly through a determination of dz/dt. In the proce- 
dure of calculating the differential ages, Simon et al. |15j 
have employed the new released Gemini Deep Deep Sur- 
vey [lal and archival data [l3l, ^M, to determine the 9 
numerical values of H{z) in the range < z < 1.8. 

The probability distribution for the parameters ai , 02 
is (sec e.g. |19| ) 



P(ai,a2)=AAe-'^'("i'"^)/2, 



(25) 



where TV is a normalization constant, 
mize the x^ function 



Hence, we mini- 



X^(ai,a2) 



E 



[Hth{ai,a2] Zj) - Hob{zi)Y 



(26) 



where Hob{zi) is the observed value of H at the red-shift 
Zi, a{zi) is the corresponding la uncertainty, and the 
summation is over the 9 observational H{zi) data points 
at red-shift Zi ^^. We adopt the priors Hq = 72 km 
s~^ Mpc~^ (mean value of the results from the Hubble 
Space Telescope key project), ri° = 0.25, fij = 0.70 and 
fi3 = 0.05 [201. Furthermore, we assume a spatially flat 
Universe with a pressureless (dust) dark matter, i.e., 71 ~ 
1. Besides, we take the theoretical expressions for the 
Hubble parameter Hth with parameters ai = r^, 02 = /3 
in the F constant case and ai ~ rrx, and a2 =72, in the 
F variable case. 

The best fit parameters aibf and a2bf are obtained 
so that Xmini^ibfiCL2bf) corrcspouds to the minimum of 
X^(ai,a2). liXmin{0'ibf,a2bf)/{N-n) < 1 the fit is good 
and the data are consistent with the considered model 
H{z;ai,a2)- Here, N is the range of data set used and n 
is the number of parameters [l9| . 
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FIG. 2: Variable F: roa vs. 72. 
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The variable x^ is a random variable in the sense that 
it depends on the random data set used. Its probability 
distribution is a x^ distribution for N — n degrees of free- 
dom. In our cases, this implies that 68% of the random 
data sets will give a x^ such that 



X^(ai, 02) - X^ {aibf , a2bf) < 2.3. 



(27) 



This equation defines a closed elliptical curve around aibf 
and 026/ in the bi-dimcnsional parameter space. The cor- 
responding Icr range of the parameter a^ is the range of 
fli for points contained within this elliptical curve. Simi- 
larly, it can be shown that 95.4% of the random data sets 
will give a x^ such that 



X^(ai, 02) - X^{aibf,a2bf) < 6.17. 



(28) 



Again, this last equation defines an elliptical curve in 
parameter space and the corresponding 2a range of the 
each parameter a^ is the range of a^ for points contained 
within this elliptical curve. 

The local minimum in the case of F constant is x^ = 
9.02445 for r^o = 0.0230767 and /3 = 0.302457. Hence, it 
follows from cq. ^ with 71 = 1 that 72 = 0.0800689 
and a — 0.899181. In the F variable case, the minimum 
is x^ ^ 9.03738 for r^ = 0.0851298 and 72 == 0.0328392. 

In Figs. 1 and 2 wc have plotted the confidence re- 
gions, that is, the sections of the elliptical curves ex- 
plained above, in the (3 vs. Too plane for a constant F 
and in the Too vs. 72 plane for a variable F, respectively. 
The points inside the inner ellipses or between both el- 
lipses identify the true values of parameters with 68.3% 
or 95.4% of probability, respectively. Also, in both fig- 
ures, the best fit value for each model is represented by 
a dot and the shady zone means the forbidden negative 
values. 



FIG. 3: Top: density parameters as functions of the red- 
shift; bottom: ratio between dark matter and dark energy as 
function of the red-shift. Straight lines - constant F; dashed 
lines - variable F. 



IV. COSMOLOGICAL SOLUTIONS 

In this section we search for cosmological solutions of 
the model proposed in previous sections. First we analyze 
the density parameters which are plotted as functions of 
the red-shift in the top of Fig. 3, the straight lines corre- 
sponding to the constant F{r), whereas the dashed lines 
to the variable F{r). One can infer from this figure that 
the energy transfer from the dark energy to dark mat- 
ter is more pronounced for the variable case, since for 
this case the growth of the density parameter of the dark 
matter and the corresponding decay of the dark energy 
with the red-shift are more pronounced than those for a 
constant F{r). This behavior can also be verified from 
the bottom of Fig. 3 which represents the evolution of 
the ratio of the two energy densities r = pi/p2 with the 
red-shift. This last figure also shows that in the future, 
i.e., for negative values of the red-shift, there is no differ- 
ence between the two cases, since both tend to a small 
values, indicating a predominance of the dark energy in 
the future. 

In Fig. 4 the deceleration parameter q = 1/2 + 3p/2p 
is plotted for the two cases. The present values of the 
deceleration parameter ^(O) and the value for the red- 
shift zt where the transition from a decelerated to an 
accelerated regime occur are : (i) q{0) ~ —0.47 and zt « 
0.72 for the constant F and (ii) q{0) ~ —0.52 and zt « 
0.68 for the variable F. These values are of the same 
order of magnitude of the experimental values: ^(O) = 




FIG. 4: Deceleration parameter as function of the red-shift. 
Straight line - constant F; dashed line - variable F. 
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FIG. 6: Ho as function of the red-shift z. 
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FIG. 5: lUi and W2 as functions of the red-shift z. Top: wl; 
bottom: w^. 



-0.74 ± 0.18 (see UM) and zt = 0.46 ± 0.13 (see H). 

The equations of state for dark matter and dark energy 
can be written in terms of the effective parameters wl = 



11 



1 and 



72 



1 as 



pi = wipi, 



VI =W2P2, 



(29) 



respectively. In Fig. 5 we have plotted the effective 
parameters wl and Wj as functions of the red-shift z. 
For the variable case, one can observe that in the red- 
shift range —l<z<2 the effective parameters as- 
sume the values between: —0.89 ^ wf ^ 7.2 x lO"'* and 
—0.89 ^ w| ^ —0.97. Hence, in this model the effective 
equation of state of the dark matter constituent behaves 
practically as a pressurcless (dust) fluid for z > and 
as quintessence for z < 0. The dark energy in the same 
range behaves always as quintessence. 



In Fig. 6 it is represented the difference between the 
apparent magnitude m and the absolute magnitude M 
of a source, denoted by /xq and whose expression is 



fio = m — M = 5 log I (1 



dz' 



+ 25, (30) 



with the quantity between braces representing the lumi- 
nosity distance in Mpc. The circles in Fig. 6 are exper- 
imental data for super-novae of type la taken from the 
work [2^. It is interesting to note that all /io - curves 
corresponding to the two models analyzed in the present 
work and the ACDM model practically coincide, the only 
small difference between them is for higher values of the 
red-shift. Furthermore, there is a good fitting of these 
curves with the experimental data. 



V. CONCLUSIONS 

To sum up: in this work a cosmological model for 
the present Universe was proposed whose constituents 
were the non-interacting baryonic matter and interact- 
ing dark components. The interaction between the dark 
energy and dark matter was related with their effective 
barotropic indexes, which were considered as functions 
of the ratio between their energy densities. Furthermore, 
the parameters of the functions were adjusted by consid- 
ering best fits to Hubble function data. Two functions 
were analyzed, namely, a constant and a variable one. 
It was shown: (i) for the proposed functions the ratio 
r is asymptotically stable; (ii) the energy transfer from 
the dark energy to the dark matter is more efficient for 
the variable case; (iii) for both functions the dark en- 
ergy density predominates in the future; (iv) for both 
functions the present values of the deceleration param- 
eter and the values of the transition red-shift - where 
the regime changes from a decelerated to an accelerated 
- are of the same order as the experimental values; (v) 
for the variable function, the effective ratio between the 
pressures and energy densities for the dark components 
indicates that the dark energy behaves as quintessence 
whereas the dark matter as practically a pressurcless fluid 
(dust) when z > 0; (vi) the behavior of the parameter ^o 



- related with the luminosity distance - with the red- 
shift does not have a sensible difference for the two func- 
tions and practically coincide with the one of the ACDM 
model, indicating a good fit with the experimental values. 
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